Liquid-vapor transition of systems with mean field universality class 
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We have considered a system where the interaction, v{r) = vjs{r) + £,'^VMF{r), is given as a 
linear combination of two potentials, each of which being characterized with a well-defined critical 
behavior; for uis(r) we have chosen the potential of the restricted primitive model which is known 
to belong to the Ising 3D (IS) universality class, while for VMF{r) we have considered a long-range 
interaction in the Kac- limit, displaying mean field (MF) behavior. We study the performance of two 
theoretical approaches and of computer simulations in the critical region for this particular system 
and give a detailed comparison between theories and simulation of the critical region and the location 
of the critical point. Both, theory and simulation give evidence that the system belongs to the MF 
universality class for any positive value of ^ and that it shows only non-classical behavior for ^ = 0. 
While in this limiting case theoretical approaches are known to fail, we find good agreement for the 
critical properties between the theoretical approaches and the simulations for larger than 0.05. 
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I. INTRODUCTION 

The exact determination of the critical properties of 
a fluid (i.e., in terms of the location of its critical point 
and of its universality class) represents a formidable chal- 
lenge both to theoretical approaches as well as to com- 
puter simulations. Although integral-equation theories 
are quite successful in predicting the thermodynamic and 
structural properties of a variety of simple and more com- 
plex systems over a wide domain in temperature-density 
spaced they meet with variable success in the critical re- 
gion due to an unsatisfactory treatment of long wave- 
length fluctuations which are of particular relevance in 
this region^. To give an example, for the modified hy- 
pernetted chain equation (MHNC), which is amongst the 
most accurate liquid state theories^, the boundary of the 
density region at which no solution is found does not 
correspond to the spinodal and there is no divergence of 
the compressibility. The mean spherical approximation 
(MSA), on the other hand, satisfies scaling laws in the 
critical region, though with peculiar critical exponents 
(mean spherical exponents^.) , but, similar to MHNC, fails 
to give a proper treatment of the first order gas-liquid 
(GL) transition^. For computer simulations, on the other 
hand, suitable techniques have been developed that allow 
the determination of the critical point via data extrap- 
olation from simulations performed in finite simulation 
volumesSiLS. 

The purpose of the present paper is to examine the 
performance of theoretical approaches and of computer 



simulations in the critical region for a system where the 
potential v{r) is a linear combination of two interactions, 
exhibiting different, but well-established critical behav- 
ior, i.e., Ising 3D (IS) and mean field (MF): 



(1) 



The dimensionless parameter indicates the relative im- 
portance of the two contributions. For reasons given be- 
low, we shall use a Coulombic interaction for vis (r) and a 
Kac potential^' for wmf('') which we define as the limit, 
at fixed volume V, 



VMF{r) 



lim a 
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r/(T, a length scale, q has 

the dimension of a charge). When the dimensionless pa- 
rameter a* — aa tends to zero the strength of the po- 
tential decreases and its range increases to infinity. It 
is crucial to observe that the limit a ^ in Eq. 
must be taken before taking the thermodynamic limit, i. 
e. at fixed volume. Although one might expect that the 
critical behavior of the system is the result of a competi- 
tion between the Ising 3D and the MF behavior, leading, 
for example, to a cross over between the two universal- 
ity classes, it turns out that the criticality of the system 
is described by a MF behavior except at = where 
the system behaves Ising 3D like. This is unambiguously 
confirmed both by simulations and theory (cf. Sees. Ill 
and IV). 
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To overcome the problem of liquid state theories men- 
tioned above, two advanced liquid state theories have 
been proposed in recent years that are able to pro- 
vide accurate data in the critical region: the hierar- 
chical reference theory (HRT) of Parola and Reattc"^ 
and the self-consistent Ornstein-Zernike approximation 
(SCOZA) originally proposed by H0ye and Stellii*i^. 
Compared to standard liquid state theories they are 
more successful in describing the critical region, provid- 
ing, in particular, non-classical (i.e., non-MF) critical ex- 
ponents. The HRT approach is a successful mergeri^i^ 
of liquid state theory with rcnormalization group ideas. 
In the simplest form of SCOZA-^-^'ij^, (which will be 
used here) the Ornstein-Zernike (OZ) equation is supple- 
mented by an approximate, MSA-type closure relation 
for the direct correlation function containing a density- 
and temperature-dependent function which is determined 
by imposing self-consistency between the compressibil- 
ity and internal energy routes to the thermodynamics. 
SCOZA predicts non-MF critical exponentsi^; they are 
different when approaching the critical temperature from 
above or below and, in both cases, satisfy hyperscalingi^. 
We note that the formalism of SCOZA is not applicable 
for the Kac potential (i.e., for a = 0) itself; instead, we 
rather have to consider a potential with a small, but fi- 
nite parameter a. On the other hand, this restriction 
has also the attractive feature that we can study in de- 
tail how the critical behavior varies as a tends to zero; 
indeed we can identify a cross-over behavior between a 
non-classical and a MF behavior. 

In addition we use a recent improved mean field the- 
ory based on a loop expansion of the free energy of the 
RPMiLi^ to describe the thermodynamic properties of 
our system. 

Finally, we have used grand canonical Monte Carlo 
(GCMC) simulations performed on a hyperspherei2iS£ to- 
gether with histogram reweighting techniquesSi and finite 
size scaling analysis to locate the critical point and the 
near critical coexistence curv o^i^^ . 

With respect to the contribution to v{r) with Ising 3D 
criticality we resorted to the restricted primitive model 
(RPM). This choice has three advantages: first, it is 
now well-established by simulations that the RPM be- 
longs to the Ising 3D universality class2i24. In addition, 
the critical parameters (i.e., temperature Tc and density 
Pc) and the near critical coexistence curve are accurately 
known. Second, the reduced, dimensionless, critical tem- 
perature and density of the RPM « 0.0489 - 0.0492, 
Pc « 0.076 - O.08O2i212^ (for the definition of the re- 
duced units see Sec. II) differ notably from those of the 
Kac model (lim^2^o T'c/C^ ~ 1.13052, pc w 0.27) which 
would not be the case by taking for wis('') a short range 
Yukawa potential. Third, the improved MF theory out- 
lined above can be formulated with closed expressions for 
a size-symmetric (and possibly charge-asymmetric) sys- 
tem of charged hard spheres, including thus the RPM—. 
These attractive features are contrasted by one serious 
disadvantage: theoretical approaches that have been pro- 



posed up to date in the literature are not able to provide 
results for Tc and pc for the RPM that are in reasonable 
agreement with the predictions of simulationi2i2Si2L2&. 
Nevertheless, the critical exponents we obtain from our 
SCOZA investigation for the RPM show a cross-over be- 
havior to MF. 

The paper is organized as follows: after presenting our 
model (Sec. II), we give a brief introduction to the theo- 
retical concepts that we apply to study the critical behav- 
ior of our system: an improved MF theory and SCOZA 
(Sec. Ill); the section is closed with a comparison be- 
tween the data produced by these concepts. In Section 
IV we give details about the simulation techniques that 
we apply and discuss the problems hereby encountered. 
In the subsequent section (Sec. V) we make a comparison 
between the data obtained in the theoretical approaches 
and in the simulations. The paper is closed with con- 
cluding remarks. An appendix provides further technical 
details about our reweighting procedure. 



II. MODEL AND REDUCED UNITS 

Using the RPM for vis{r), our potential u(r) is actually 
the interaction of a binary system; thus 

{oo r < <T 
, e2 / N ^ . (3) 
1- VMFir) r > a ^ ' 
r 

where VMF{r) is given in (j21). 

Numerical results will be expressed in reduced units 
(taking the hard-sphere (HS) diameter a as unit of 
length): reduced temperature T* = 1//3*, where (3* = 
q^/k-oTa {q — q+ — — q_ charge of the spheres, fee 
Boltzmann constant, T temperature), reduced configu- 
rational chemical potential p* = p/k^T, reduced volume 
V* = V/a^ and reduced density p* = N/V*{N the num- 
ber of particles). However, for notational convenience we 
will drop the stars throughout the paper. 



III. THEORY 

A. Improved mean field theory (MF2L) 

The Kac potential energy Ep(C) of a grand canonical 
(GC) configuration C, assuming, for instance, periodic 
boundary conditions, is 

^p(^) = ^ E vMFik)mpi-k) ; (4) 

where VMF{k) and p(fc) are the Fourier transforms of 
'^mf(^) and the microscopic density, respectively. For 
a fixed volume vuFik) = ^{k/a) vanishes for all non 
zero (fc 7^ 0) wavevectors of Fourier space in the limit 
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a ^ and umf(O) = 3>(0). Thus 

1 

W 

-2TTq^Vp' 



(5) 



The GC partition function of the model Eq. Q is then 
obtained asi^ 

/•OO 

Sy (/?, I/) ~ / dp <T^ exp [-V {-vp+ 
Jo 

+PfRPM{p)-2nl3q^eP^)] , (6) 

where we have replaced the sum over the number of par- 
ticles by an integral over the density which is valid for 
large systems and has no consequences for the thermo- 
dynamic hmit (here (3 = I/^bJ") . In Eq. /rpm(p, P) 
is the Hclmholtz free energy per unit volume of the RPM 
fluid at density p and v = (3 p.. Taking now the limit 
V ^ OO, one gets 

lim - In Sy (/?,;/)= min£(/3, 1/, p) , 

C{P, V, p) - -up + pfnpuip, P) - 27rPq^eP^ ■ (7) 

The function C{f3,u,p) plays the role of a Landau 
functio n--'" . Denoting /3^^^ the inverse critical temper- 
ature of the RPM, we note that /3/rpm(p, f3) is an analyt- 
ical and strictly convex function of p for /3 < Pf'^^ and 
that the critical point of our model occurs for (3 < /?^^^ . 
From the analysis of the Kac model it can thus be in- 
ferred that the composite system Eq. will have MF 
behavior^'-'^®. 

Equation Q can be taken as the starting point for 
an improved mean field (MF) theory. To this end an 
approximate expression for /rpm(/3,p) is inserted in the 
Landau function C{f3, u, p) of Eq. Q and C{(3, u, p) min- 
imized with respect to p for each inverse temperature (3 
and reduced chemical potential v. At low temperatures 
and arbitrary v, C{(3, v, p) has in general two minima, pg 
and pi, but there is only one value of v, i.e. i'cocx(/3), for 
which the Landau function takes equal values at the two 
minima. The corresponding minima are the densities of 
the gas and liquid phases at coexistence, respectively. 

We have used for /rpm the loop-expansion of the RPM 
free energy obtained in the field theoretical framework of 
Rcf. 17,18.1 . i.e. 



f _ f(0) I f(l) 

/RPM — /rpm + /rpm 



f(2) 

/rpm 



(8) 



where the superscripts '^^^ denote the p*"^ order contribu- 
tion to the loop expansion of /rpm- 

At the tree level one has for the RPMii 



a^/rpm 



/3/Hs(p)-pln2- ^/3gV(0) , 



(9) 



where /hs (p) denotes the excess free energy per unit vol- 
ume of the reference HS fluid for which the Carnahan- 
Starling approximation^ was used in numerical calcula- 
tions. The pin 2 contribution stems from the entropy of 



mixing and Vc{r) denotes the Coulomb potential, Vc{r) 
1 /r for r > a and regularized inside the core region. 
The one-loop contribution f3fl^pM reads 



PL 



(1) 

RPM 



(2^ 



\n{l + I3pq^vc{k)) , (10) 



and coincides with the free energy in the random phase 
approximation (RPA)i2,. 

(2) 

Finally, the two-loop contribution PfjipM^ ^^^^^ 
tractable, has expression 



(2) _ 0^ 



RPM 



^[pq'rj d^r/^HS,p(r)A2(r), (11) 



where /ifjs, pi''') denotes the usual pair distribution func- 
tion of hard spheres (HS) at number density p; A(r) the 
propagator of the free theory whose expression in k- space 
isU 



A(fc) 



Wc(fc) 



(12) 



+ f3 pq^ Vc{k) 

As discussed at length in Ref. PIT^ the loop- 
expansion is not a systematic expansion in some small 
physical parameter and thus depends explicitely on the 
regularization of the Coulomb potential Vc{r) in the 
core (r < a). In this work we adopted the MSA 
regularization^^ which ensures that, at the one-loop 
level, the pair distribution function vanishes inside the 
core. Moreover, with this specification the one-loop ap- 
proximation of /rpm coincides with the optimized RPA 
(ORPA) theory or the MSA theory if one adopts, as we 
did, the Percus-Yevick (PY) expression^ for /ihs, p{f)- 

With the MSA regularization^ 



VcMSAir) = l/r 

B f Br 



B 



{r>a), (13) 



(0 < r < cr) , (14) 



c2+a;-.T(l + 2a;)i/2 



(15) 



where x = na and = An^pq^ is the Debye wavenumber 
squared. 



B. SCOZA 

As noted in the introduction the formalism of SCOZA 
is not applicable if we choose in our interaction Vij (r) 
a to be zero. Therefore we consider a slightly modified 
potential, Vij{r), the so-called charged Yukawa model^ 

«y W - v?A^) + ^^^^(r) = +' - (16) 

where vf^ (r) stands again for the RPM while the Yukawa 
(Y) contribution reads 



00 



v^[r) 



-q e 

r 



-Q*r/fj 



r < a 
r > a 



(17) 
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In the limit a* 0, (r) reduces to vupir)- 

For this particular system the set of OZ equations, 



hij{r) = Cij{r) + '^Pk [hik Ckj] {r 



hi, k ■ 



(18) 

where '(g)' denotes the convolution operation, can be de- 
coupled: with ^ p-\- = p/2 and q+ — —q~ — q and 
introducing the correlation functions 







fM 


+ 




(19) 


c^ir) = 


1 r 


hM 


+ C+ 


-W] 


(20) 


h^ir) = 




f(r) 


-h. 


-M] 


(21) 






hM 


- c+ 


-M] 


(22) 



the set of OZ equations (|18() decomposes into two inde- 
pendently solvable one-component OZ equations for the 
Yukawa and the Coulomb parts, namely 



h^(r) = c^(r) -t- p 
h'^lr) = c'^(r) -f' p 



(r) 
(r). 



(23) 



The formahsm of SCOZAli^ is based on the MSA 
closure to the OZ equations, i.e.. 



9ij (^) ~ ^ r < a 

Cij{r) = -I3v^j{r) for r > a. 



(24) 



Introducing [r) — (r) + 1 and g''^{r) — hP{r), leads 
to the MSA closure for the decoupled OZ equations ()23|l : 



h^{r) = -1 
cY(r) = 



r < a 



(r) — (3v^ {r) r > a 



and 



/iC(r) 
cC(r) 



7' < (T 
7' > (7 



(25) 



(26) 



Above we have used for the correlation function of the 
HS reference system, c^^(r), the parameterization (for 
details see Ref. js^) 



„HS 



(r) = K,{p) 



exp[-zi{p)(r - a)] 



(27) 



For the Yukawa part the formalism of SCOZA can be 
applied in a straightforward way; introducing a yet un- 
known, state-dependent function K^{p,T), the SCOZA 
closure relation (I25II reads 



h^{r) 



-1 r < a 

(r) = c"S(r) - K^{p, T)v^{r) r > a. 



(28) 



For the Coulomb part the situation is more delicate: in 
Coulomb systems the Stillinger-Lovett sum-rulea^*^ have 



to be satisfied which would be violated when introduc- 
ing a function K'~^{p,T), similar to what we did for the 
Yukawa contribution in 1)28(1 . We therefore cannot ex- 
tend SCOZA to the Coulomb part and have to treat it 
rather within MSA, i.e., we use the closure relations H26() 
for h'~^{r) and c^(r). 

Within the SCOZA formalism a partial differential 
equation (PDE) is derived that imposes thermodynamic 
self-consistency between the energy and compressibil- 
ity routes and thus fixes the yet undetermined function 
K^{p, T). Thus we have to calculate the internal energy 
and the compressibility of the system within the frame- 
work of our theory. A straightforward analysis leads to 

M = — = 27r ^ pipj j drr^gij{r)vij{r) ^ + 

(29) 

with u^^ being the excess (over HS) internal energy and 
introducing 



27: p^ I drr^ {r)v^ (r) 



(30) 



and 
,c 



/3u^= 27rp2/3 / g'^{r)q^rdr (Vl + 2w - 1 - w) , 

Jo 47ra-^ 

w = ^/ At: pa^q'^/kBT. (31) 

The explicit expression of was derived within MSA 
by Waisman and Lebowitaffi. 

In a similar manner the reduced, dimensionless com- 
pressibility Xrcd(= P^^Txt) can be calculated for our 
system, leading to 

(Xrcd)"' = 1 - ^ [5++(0) + c+_(0)] = 1 - p5Y(0) (32) 

where the tilde denotes the Fourier transform. 

SCOZA imposes thermodynamic consistency between 
the energy and the compressibility routes via the follow- 
ing PDE 



1 



dp^ df3 Vxrcd 



(33) 



Using (|29|l - (|32|l and a few trivial manipulations one 
ends up with the following equation 



B{p,u-f-^+D{p,u^)^C{pf^, 



(34) 



i.e., a PDE which has a structure similar to the standard 
SCOZA-PDE for the one-component fluid (see e.g. Ref. 
^3). B{p,u^) is formally identical to the corresponding 
standard SCOZA-expression'^® B{p,u), and is obtained 
for the present problem by simply replacing u by u^; 
again, C{p) = p. Finally, one finds for the additional 
coefficient D{p,u^) the following relation: 



D{p,u^ 



2,,C 



(35) 
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Once the PDE H34() is solved, further thermodynamic 
quantities (notably the pressure and the chemical poten- 
tial) follow via standard thermodynamic relationsi^*^. 

While previous applications of SCOZA were rather 
dedicated to the determination of the coexistence curves 
and to a reliable location of the critical point within rea- 
sonable accuracy, we rather focus in the present contri- 
bution on a quantitative description of the critical prop- 
erties. This means that we have to approach the critical 
point very closely: introducing r = {T — Tc) /Tc, we need 
accurate and reliable data down to r ^ 10^^° or even less. 
Compared to previous SCOZA applications this requires 
a considerable increase in the numerical accuracy. We 
have realized this goal by suitable modifications of the 
code and by increasing the default numerical accuracy 
to four-fold (extended) precision. Taking these measures 
is in particular indispensable for small a- and ^^-values 
where the Coulomb and Yukawa potentials at contact dif- 
fer by many orders of magnitude. Of course this increase 
in numerical accuracy is accompanied by a considerable 
increase in computational effort. 



C. Comparison between the theoretical approaches 

In Figure ^ we show the phase diagram of the system 
for = 1. We display results obtained from the im- 
proved MF theory both on the one- and on the two-loop 
level. We show also SCOZA results, varying the index a 
from 1.8 down to 0.01, which - within numerical accu- 
racy - is already very close to the Kac- limit. We observe 
that in this limit the SCOZA data really 'converge' to 
the results obtained within the improved MF theories. 
A more detailed analysis shows that for a = 0.01 the 
SCOZA data coincide with the results obtained on the 
one-loop level, while there is a slight difference w.r.t. to 
the two-loop level. This outcome is valid for any value 
of ^. Although a formal proof may be difficult, the limit 
a — > of SCOZA will reproduce the exact MF result 
for h^{r) and c^(r) while h'~^{r) and c^(r) are solution 
of the MSA. The one-loop theory is precisely built with 
these ingredients and should coincide with SCOZA as 
a ^ 0. The MF2L theory improves h^{r) and c^(r) for 
all ^ although for ^ — > the improvement is minor (cf. 
Sec. V). Comparison with simulation data will show that 
the two-loop level indeed represents an improvement. An 
even more elaborate version of the two- loop version of the 
improved MF theory that requires the pair distribution 
function to vanish inside the core will certainly lead to 
an even better agreement with simulations. 

In Figure 121 we show results obtained for the effective 
critical exponent 7cff as a function of r, obtained from 
SCOZA for = 0.16 and for a sequence of a- values; a 
varies from 1.8 down to 0.01. 7 characterizes the diver- 
gence of the compressibility as one approaches the critical 
point from above, i.e., xt . The effective expo- 

nent 7cff is obtained by differentiating the logarithm of 
XT w.r.t. the logarithm of r along the critical isochore. 



While in the MF universality class 7 is 1, its SCOZA 
value^^ is 2. Rather far from the critical point (i.e., for 
T ~ 10^^) we observe a MF behavior; as we approach 
the critical point, a cross-over to the SCOZA behavior 
takes place. In addition, the region where the SCOZA- 
value is attained shrinks drastically as a is lowered. This 
reflects that in the Kac-limit a ^ the MF behavior be- 
comes dominant, just as we expect it to be. The curves 
presented in Figure |21 indicate that for a = we would 
observe a MF behavior over the entire r-range. 

IV. SIMULATIONS 
A. Methods 

A detailed test of the validity of the theoretical ap- 
proaches has been made by comparison with GCMC sim- 
ulations. The GCMC simulations are performe d on a 
hypersphere using the approach, detailed in Ref. |l9l20j |. 
where the system of charged spheres is viewed as a single 
component fluid of charged bihard spheres constrained 
to move at the surface of a four dimensional sphere. The 
practical implementation of this method, which is par- 
ticularly adapted and efficient for simulating Coulomb 
systems of size N < 10000, can be found in Ref. [s^ . 

Sampling of configuration space of the system is made 
according to the GC probability distribution 

p(A^,i?^,C,Ai,T,y) = (36) 
^ exp {p[fiN - t/R.PM(i^^) - UuFjN, V, Q]} 
7V!A3^S(r,t/,/i,C) 

Here S is the grand canonical partition function^- and 
indicates the positions of the N charged spheres (or 
charged bihard spheres); Urpm{R^) is the energy of the 
RPM comprising Coulomb and hard core interactions. 
The term UupiN, V,^ = -aq^^^N^ (with a = 2'k/V) is 
the MF contribution to the internal energy as explicited 
in Sec. IIIA. 

For each value of ^, we determine the critical point and, 
in its vicinity, the coexistence line of the GL transition. 
The location of the latter is obtained from the histograms 

i7(iV,u,e,Ai,T,y) = (37) 
= C j dR"" 5[u-U^pm{R'')\ p{N,R^,^,fi,T,V) 

where C is a constant of normalisation. By integration 
of H{N,u,£,, iJ,,T,V) over u, one obtains the histogram 
h{N, ^, ^,T,V) which shows, in the vicinity of the GL 
coexistence line two peaks located at = Ng and Ni , re- 
spectively, and corresponding approximately to the den- 
sities pg and pi of the gas and the liquid. 

At given T and ^, the values of Pe{L) {L characteristic 
size of the system defined as V^^^) corresponding to GL 
equilibrium are affected by finite size effects. Therefore 
the estimates of pe (L) are made for increasingly large vol- 
umes and are then extrapolated to the thermodynamic 
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limit. Several procedures are possible to define for a finite 
volume V the states of GL coexistence: First one can de- 
termine He{L) such that the two peaks in h{N, ^, /i, T, V) 
have equal height or, preferably, equal area. This pro- 
cedure, however, does not allow, at fixed V, a very pre- 
cise location of the critical point as in its vicinity the two 
peaks of h{N, ^, fi, T, V) tend to merge. A second method 
to locate, at fixed V, phase equilibrium as well as critical 
temperatures and densities has been proposed by Bruce 
and WildingS. It relies on the use of histograms of the 
variable M ^ p - s{u - aq^(^N'^)/V 

p{M,^,p,T,V)^ I duY,5{M-p+ (38) 

+s {u - aq^eN^)/V) H{N, u, ^, fi, T, V) . 

For an appropriate choice of the parameter s, M is such 
that at GL equilibrium and near the critical point the dis- 
tribution p{M, ^, /i, T, V) satisfies the symmetry relation 

p{M- (M), t M, T,V)^p i-M + (M) , Ai, T, V) 

(39) 

expected for the order parameter that characterizes 
phase transitions of the Ising 3D or the MF universality 
class. It is remarkable that, quantitatively, the histogram 
h(N, ^, /i, T, V) corresponding to the symmetrized distri- 
bution p{M,^, fi,T,V) has two peaks of equal area. A 
third possibility is to apply an unbiased finite size scaling 
method presented by Fisher and co-worker32i2&. 

1 /2 

If one defines x — {SM) / (SM'^) with 6M ~ 
M. — {M), the critical temperature Tc{L) is obtained by 
determining the values /x, T, and s which enable to fit 
p{x,^, fi,T,V) = p{x) to the critical distribution Pc{x) 
of the order parameter of a system of known universal- 
ity class. For example Pcix) can be the distribution of 
magnetization of the Ising 3D model or that of the order 
parameter of a lattice spin model of the mean field uni- 
versality class. For the Ising 3D model Pc{x) is known 
accurately from MC simulations''^ and exactly for the 
MF universality class (see below). The values of Tc{L) 
obtained for increasingly large volumes V can then be 
extrapolated to obtain Tc(oo). 

From the histograms p{x) obtained at fixed volume 
the moments {6A4'^') and (JA^'*) can be calculated 
and from these the fourth order cumulant, Ql(/3) = 
(SM'^) I {SM'^f, which in the vicinity of Tc{oo) = 1/(3^ 
has the formffi 

Qdfi) = Q* + ai (/3 - /?,) + a2 (/? - Pcf i'"^ + 

+03 {13 - Pcf L^"- + ... + biLy^ + ... . (40) 

Q* and the exponents yr and yi are characteristic of the 
universality class of the system and can be determined, 
together with /3c, by a fit of the Ql{I3) obtained for dif- 
ferent system sizes. 

For the MF universality class p{x) = Aexp {—ax'^) 

with A = [r(3/4)]22l/4/^3/2^ ^ ^ ljp(3/4)]4/^2^ Q* ^ 



2[r(3/4)]V7r2 = 0.456947, = d/2 = 1.5 and = 
— (i/2 = —1.5, respectively, where d is the dimensionality 
of the system^ii^. As is well known the MF exponents 
do not satisfy hyperscaling. 



B. Results 

Simulations have been performed for the values = 
100, 11.11, 1.0, 0.1, 0.05, and 0.02, and volumes V = 
1000, 4000, 8000, and 16000. For each pair of f and V 
of the order of 5 — 10 x 10^ MC trial moves were sampled. 
At volume V = 16000, even for such a large number of 
configurations the statistical error remained large so that 
the determination of phase equilibrium and the location 
of the critical point was based solely on the smaller vol- 
umes. Further, a reweighting procedure allowed to ex- 
tend the results obtained at the ^-values listed above to 
neighboring values of ^ (see Appendix). 

An approximate location of the GL transition inferred 
from the appearance of two peaks in the histograms 
h{N,i^, p,,T,V) indicates that {Tc, Pc) varies from ap- 
proximately (113, 0.25) at = 100 to (0.075, 0.11) at 

= 0.02. From these values we can conclude that for 
the range of ^ values considered the critical temperature 
of our model is markedly higher than the critical temper- 
ature of the RPM (^2 ^ 0), i.e., - 0.0492^. 

It follows from the analyticity of the free energy (Sec. 
IIIA) that the critical behavior of our model can be an- 
alyzed according to the classical scheme (see e.g. Ref. 
[43j)- In particular, the critical exponents derived from 
the fit of Ql iP) should be MF-like and in the procedure 
of Bruce and Wilding£, p{x) has to be adjusted to the 
distribution Pc{x) of a model with MF criticality. 

In Table B we summarize the critical temperatures 
estimated from the fit of Ql{P) calculated for the vol- 
umes V = 1000, 4000, and 8000 at the six values of 
considered in the simulations. These temperatures 
are in good agreement with values for Tc{L) obtained at 
F=8000 by fitting p{x) to the distribution Pc{x) corre- 
sponding to the MF universality class. In the table Pc{L) 
denotes the average density of the thermodynamic states 
at Pc{,L) and Tc(L) corresponding to these fits. 

The considered values of were too different to allow 
reweighting in the range of parameters > 0.1. How- 
ever, it was possible to combine histograms obtained at 
^^=0.05 and 0.1, on one hand, and those at — 0.02 and 
0.05, on the other. Furthermore, the results at = 0.10 
could be reweighted up to 0.13, and those at — 0.02 
to 0.01. However, combination of the H{N, u, ^, p, T, V) 
was possible only at volume V — 1000 due to insufficient 
overlap of the energies at the larger volumes. 

Figures 13 and ^ show the variation of Tc{L) and Pc{L) 
as a function of between ^=0.01 and 0.13 for V = 1000. 
GL coexistence curves are presented in Figs. |5l - |H1 for 
C^=0.02, 0.05, 1.0, and 11.11 and the three volumes 
V =1000, 4000, and 8000. The error in Tc{L) is of the 
order of 0.3 — 0.1% depending on the amount of sam- 
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pling of the histograms h{N,£^, iJ,,T,V). Tc{L) decreases 
weakly (~ 1%) when V varies from 1000 to 8000 and 
differs from by about 1%. From these results it is 
concluded that in the range 0.01< <0.13 the values 
of Tc{L) obtained by reweighting for V =1000 provide 
a rehable estimate of or Tc{L) at V^=8000 with an 
uncertainty of the order of ~ 1% (cf. Table II). 

The variation of the critical density with volume is 
small (cf. Table II), of the order of the statistical error 
on h{N, ^, ficiL), Tc{L), V), estimated to be ~ 0.5 - 1.0% 
depending on sampling in the simulations. A second es- 
timate of pc, denoted pi{L), is based on the rule of rec- 
tilinear diameter, i.e., p^(L) = [pg{T, L) + pi{T, L)]/2 
valid at a MF critical point for T close to Tc- The values 
of pf{L) given in Table II for the different values of ^ 
and V , were obtained by retaining only the coexistence 
densities pg{T,L) and pi{T,L) which were closest to the 
critical point at Tc{L) and Pc{L). The coexistence densi- 
ties Pg{T, L) and p;(T, L) correspond to the two maxima 
of h{N,£_, p,T,V) whose associated distribution p{x) is 
symmetric (see above). The uncertainty of the location 
of the maxima of h(N,^, p,T,V) increases when T ap- 
proaches Tc{L) so that the error on p^iL) is of the order 
of ~ 2%. The result for Pc{L) is somewhat lower than 
Pc{L) though compatible within the combined error esti- 
mates. 

As mentioned previously in the literature^, agreement 
ofp{x) with aiPc{x) characteristic for a universality class, 
is necessary but not sufficient to make unambiguous con- 
clusions about the universality class of a model system 
when studied by simulation. This uncertainty is illus- 
trated for ^2 = 0.02 and V = 8000 in Fig. which shows 
excellent fits of both the Pc{x)'s correponding to the Ising 
3D and MF universality classes. 

Although for the present model the MF universality 
class is implied by the analyticity in T and p of the free 
energy F{£^, T, V) in the vicinity of the critical points de- 
termined for > 0.01, this critical behavior is also con- 
firmed by the fit of the Ql (/?) obtained for the different 
values of V. To realize the fits, Ql{P) was calculated, for 
each volume, along the coexistence line by the procedure 
outlined above, and, for T > Tc{L), for values of p and s 
such that p{x) remained symmetric with respect to x = 0. 
For each volume, 8-12 values of Ql{P) were determined 
in the neighborhood of Tc(L) [i. e. « ±3% from Tc[L)\. 
The fit is obtained by minimizing the square deviation 
between the simulation results at the three volumes V — 
1000, 4000, and 8000 and Eq. ^ retaining in the ex- 
pression the eight parameters shown explicitly. In Table 
III the values of Q*, /?^, yr, and yi are summarized for 

= 0.02, 0.05 and 0.1. The precision of the fits can be 
appreciated from Figs.llOlandllll Taking into account an 
uncertainty of 1 — 2% on Ql{P, L), fits compatible with 
this error were obtained for a wide range of values of Q* 
and yi, determined, respectively, with a precision of ~ 10 
and ~ 30%. In contrast, the precision of = l/T^ is 
~ 1% and of yr ~ 5%. The value of yr is in agreement 
with the MF value and, despite large uncertainties, also 



the values of Q* and yi. 



V. COMPARISON: THEORIES VS. COMPUTER 
SIMULATIONS 

While the results obtained of the two theories and the 
simulation data have already been discussed in the re- 
spective subsections we focus in the following on a direct 
comparison between the theoretical and the simulation 
data. 

To this end we go back to Table □ where we have col- 
lected critical data (i.e., Tc and pc) for different values 
of as obtained from the GCMC simulations, the im- 
proved MF theory, and SCOZA; in the latter case we 
have chosen a — 0.01, i.e., a value that has turned out to 
be sufficiently small to make the potential {r) Kac-like. 
From these data we can conclude that both theories are 
able to reproduce the critical temperatures and critical 
densities for large and intermediate values of > 1 very 
accurately. In the range of values 0.06 — 1 the MF2L 
approximation remains quantitative while SCOZA dif- 
fers from the simulation results by « 10% for the critical 
temperature and ~ 30% for the critical density. Below 

~ 0.06 the RPM contribution to the potential becomes 
dominant and we find, not surprisingly, that the theoret- 
ical data start to differ from the simulation results; in 
both theories, critical temperature and density tend to 
the MSA values of the RPM for = q whose failure 
to reproduce the simulation data for the RPM is well- 
documentedii. As pointed out above SCOZA identifies 
with the one-loop MF theory for all as a ^ 0. 

The coexistence curves in the vicinity of the critical 
point obtained from theory and simulation are compared 
in Figs. IS]- IHl for selected values of <^^. The convergence 
of the coexistence curves with system volume can be ap- 
preciated from these figures; the results obtained for vol- 
ume V = 8000 appear to be close to the infinite sys- 
tem limit. The MF2L approximation for the coexistence 
curve agrees well with simulation data for > 0.1 but 
detoriates when approaching the RPM limit where both 
the critical temperature and density are overestimated. 
One should note, however, that even at = 0.02 the 
critical temperature differs from the MC result by only 
10%. The MF2L approximation is a definite improve- 
ment over either the one-loop approximation or SCOZA 
and could still be improved by requiring the pair distri- 
bution function to vanish inside the core at the two-loop 
level. 

As noted above, SCOZA provides data, that are in 
a less satisfactory agreement with simulation data than 
those of the improved MF theory. Despite this deficiency, 
SCOZA was able to provide useful information about the 
critical behavior of the system: having non-classical crit- 
ical exponents we have been able to trace the cross-over 
from a MF to a non-MF behavior, as we approach the 
Kac-limit (see Figure 2), verifying thus that the system 
displays MF criticality for a = 0. 
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VI. CONCLUSION AND OUTLOOK 

We have investigated the critical behavior of a system 
where the interaction is a Unear combination of two po- 
tentials, each of them having a well-defined critical be- 
havior, viz. Ising 3D or MF. Taking the RPM for the 
first contribution and a Kac-like potential for the latter 
one, we could vary via a parameter the weight of the 
respective contributions. Our investigations have been 
carried out with GCMC simulations and two theoretical 
approaches, i.e., SCOZA and the improved MF theory. 
The theoretical argument of Sec. Ill shows that MF be- 
havior at the critical point is expected for all ^ 0. In- 
deed this conclusion can unambiguously be drawn both 
from simulation and SCOZA data. 

We have given detailed account how the simulations 
and the theoretical approaches perform in the critical re- 
gion. Wc could demonstrate that both simulations and 
theory have reached a level of sophistication and accu- 
racy that allows a quantitative description of the critical 
behavior of simple liquids in terms of the location of the 
critical point as well as of the universality class. To be 
more specific, the location of the GL critical point (based 
on the MF universality class) and the coexistence lines 
are obtained from simulation within a numerical accu- 
racy of ^ 1 — 2% with a slight variation with the volume 
V. For our model we find excellent agreement between 
the theoretical approaches and computer simulations for 
^^-values down to ^ 1; the improved MF theory, MF2L, 
performs in an excellent way even down to ~ 0.05. 
Only for very small values of i.e., close to the pure 
RPM, simulation data and theroretical results start to 
deviate. Nevertheless, SCOZA predicts the correct criti- 
cal behavior, i.e., MF. 

We highlight the ambiguity arising when determining 
critical temperatures by fitting universal critical point or- 
der distribution functions by showing that for the present 
system almost equally good fits are obtained for both the 
universal Ising 3D and MF distributions. 

Finally, by referring to the theoretical developments 
presented in Sec. IIIA we can remark that the critical 
temperature is obtained by requiring that the second 
derivative of £{(3,1/, p) with respect of the density van- 
ishes which yields 



(41) 



where Xt denotes the isothermal compressibility of the 
RPM. Arguing that for small S,"^, xt ~ \Tc - Tf^^[i 
where 7 is the critical exponent of the compressibility of 
the RPM (or, more generally, that of the system with 
potential t;is) it is concluded that 



This relation provides a means to evaluate the expo- 
nent 7. By simulation this turns out to be practical only 
for large systems {V > 20000) which are necessary to ac- 
curately sample the density distribution function in the 
very low density gas phase of the RPM near the critical 
point. 
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APPENDIX: REWEIGHTING 

Reweighting starts with the histograms 

Hh{N,u, pk, Pk,^k) obtained with Afk entries at 

fixed volume V, chemical potential pk, and inverse 
temperature pk (« = 27r/F) 



Hk{N,u,pk,0k,^k) = 



J^k 



I 



S(/3fc,/ifc,&) 

A4 



= 0.{N,u)- 



Pk tikN + i3k ilaN^'-ffkU 



'^ifik,IJ'k,^k) 



(A.1) 



With n histograms (k = 1, ...n), an estimate of the 
density of states fl[N, u), up to a multiplicative constant, 
is given by 



n{N, u) 



Efe Hk{N,u,Pk,Pk,^k) 



hRPM 



2/7 



(42) 



^ e^(^^'^^'€0 Efc Hk{N,u,h,Pk,ik) 

Y^^Mk e^'' l^l'^ + l^'' 9^ ilaN^-f3kU-f{0k,n,M 

with e--^('5'= = 'B{l3k,Pk,^k) and e-f^'^'"^""^''^ = 

gF(/3i,^i,6)-F(/3fc,/ife,£fc). 

From the knowledge of fl{N,u) any histogram 
Hj {N, u,^j, pj , 13 j , V) can then be calculated up to a mul- 
tiplicative constant, according to 
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Using the initial guess /o(/3fc, Mfc? Cfc) for f{f3k,tJ'k,^k), a first estimate of Q,{N,u) is obtained from 

Efc Hk{/3k,fik,£.k,N,u) 



n 



-/3™ 9" ei)aA'2-(/3fc -/3™)"-/o(/3fc,Mfc,Cfc) 



where the factor Bm{N^u) = e P'^t^^^N p^q i^aN p^u jg introduced to prevent too large exponents in the sum 
Y,k in the denominator with = {Y.k=i,..,nf^k)/n, fi^ = {Y.k=i,..,nf^kfJ-k)/{n (3^^), and ^„ = (I]fc=i,..,„ 
A new estimate of f{PkTlJ-kT£,k) is then obtained from 

/ duY,H{N,u,Pm,^im,U) e('^'=^'=-^"'"")^ + ('''=«'~'''"«'«")"^'-(^'=-''")" = e^^(''-'''='«'-) (A.3) 

•' N 

where / du is a discrete sum over the energy values u effectively reached in the n MC runs. The values of fi{Pk, fJ'k, £,k) 
lead to a new estimate of H. Iteration is continued until for each k we find \ fi{Pk, IJ'k, £,k) — fi-i{Pk, f^k, £,k)\ < 10"''. Af- 
ter convergence H is normalized to 1. In view of the factor B„i{N, u) which relates H to i^{N, u), H{N, u, Pm, IJ'm, £,m) 
is the histogram, normalized to 1, corresponding to the average state /3„i, fJ-m, and ^m- 
From H one calculates the histograms HJl'^'^ {N,u, fik, Pk,^k) of states (3k, fJ'k, and £,k by 

Hr{N,u,^ik,Pk,^k) = H{N,u,Pr^,^i^,U) e('5'='"=-'5"'"")^ + (^^«'-'5™«'«™)"^'-('5'=-''")". (A.4) 

I 

Comparison of the histograms calculated by MC with valuable test of the convergence of the MC simulations, 
those obtained by the reweighting procedure provides a 



Electronic address: 
Electronic address: 
Electronic address: 
Electronic address: 
Electronic address: 
Electronic address: 
C. Caccamo, Phys. 



pauschenwein@cmt .tuwie nTac.atl 

jean-michel.caillol®th.u-psud.t'r| 

dominique.levesque@th.u-psud?fr] 

jean-jacques.weis@th.u-psud.fr 

elisabeth.schoeIl-paschin ger@univie.ac.at| 

gkahl@tph.tuwien.ac.at| 

Rep. 274, 1 (1996). 



A. Parola and L. Reatto, Adv. Phys. 44, 211 (1995). 
A. Parola and L. Reatto, Phys. Rev. Lett. 53, 2417 (1984); 
Phys. Rev. A 31, 3309 (1985); for a more recent overview 
see also A. Reiner and G. Kahl, Phys. Rev. E 65, 046701-1 
(2002); A. Reiner and G. Kahl, J. Chem. Phys. 117, 4925 
(2002); D. Pini, M. Tau, A. Parola, and L. Reatto, Phys. 
Rev. E 67, 046116-1 (2003). 

J.-P. Hansen and I. McDonald, Theory of Simple Liquids, 
(Cambridge Univ. Press, Cambridge, 2006), 3rd edition. 
G. SteU, Phys. Rev. 184, 135 (1969). 

N.B. Wilding and A.D. Bruce, J. Phys.: Condens. Matter 
4, 3087 (1992). 

Y.C. Kim and M.E. Fisher, J. Phys. Chem. B 108, 6750 
(2004). 

Y.C. Kim, M.E. Fisher, and E. Luijten, Phys. Rev. Lett. 

91, 065701 (2003). 

M. Kac, Phys. Fluids 2, 8 (1959). 

J.L. Lebowitz and O. Penrose, J. Math. Phys. 7, 98 (1966). 
J.S. H0ye and G. SteU, J. Chem. Phys. 67, 439 (1977). 
J.S. H0ye and G. SteU, Mol. Phys. 52, 1071 (1984). 



1^ D. Pini, G. SteU, and N. B. Wilding, Mol. Phys. 95, 483 
(1998). 

^'^ G. Kahl, E. SchoU-Paschinger, and A. Lang, Monatsh. 
Chem. 132, 1413 (2001). 

E. SchoU-Paschinger and G. Kahl, Europhys. Lett. 63, 538 
(2003). 

" J.S. H0ye, D. Pini, and G. SteU, Physica A 279, 213 (2000). 
J.M. Caillol, Mol. Phys. 103, 1271 (2005); 
arXiv:cond-mat /0409455 

J.-M. Caillol, Mol. Phys. 104, 1931 (2006); 
arXiv:cond-mat /0602205 

J.M. Caillol and D. Levesque, J. Chem. Phys. 94, 597 
(1992). 

^° J.M. Caillol, J. Chem. Phys. 99, 8953 (1993). 

2^ R.H. Swendsen, Physica A 194, 53 (1993). 

22 N.B. Wilding, Phys. Rev. E 52, 602 (1995). 

2^ J.M. Caillol, D. Levesque, and J.-J. Weis, J. Chem. Phys. 

116, 10794 (2002). 
2"' E. Luijten, M.E. Fisher, and A.Z. Panagiotopoulos Phys. 

Rev. Lett. 88, 185701 (2002). 
25 A. Panagiotopulos, J. Chem. Phys. 116, 3007 (2002). 

G. SteU, J. Stat. Phys. 78, 197 (1995). 
2^ Y. Levin and M.E. Fisher, Physica A 225, 164 (1996). 
28 O.V. Patsahan, Condens. Matter Phys. 7, 35 (2004). 
2® N. Goldenfeld, Lectures on Phase Transitions and the 

Renormalization Group, (Addison- Wesley, New- York, 

1992). 



10 



N.F. Caxnahan and Starling, J. Chem. Phys. 51, 635 

(1969). 

T. Kristof, D. Boda, J. Liszi, D. Henderson, and E. Carl- 
son, Mol. Phys. 101, 1611 (2003). 

E. SchoU-Paschinger and G. Kahl, J. Chem. Phys. 118, 
7414 (2003). 

E. Waisman, Mol. Phys. 25, 45 (1973). 

F. H. StilUnger and R. Lovett, J. Chem. Phys. 48, 3858 
(1968). 

E. Waisman and J. L. Lebowitz, J. Chem. Phys. 56, 3093 
(1972). 

E. Scholl-Paschinger, Ph.D. thesis, Institut fiir Theoretis- 
che Physik, TU Wien, Wiedner Hauptstr. 8-10, A-1040 

Wien, Austria, 2002. 

J.M. Caillol, D. Levesque, and J.-J. Weis, J. Chem. Phys. 
107, 1565 (1997). 

Y.C. Kim and M.E. Fisher, Comput. Phys. Commun. 169, 
295 (2005). 

M.M. Tsypin and H.W.J. Bloto, Phys. Rev. E 62, 73 
(2000). 

H.W.J. Blote, E. Luijten, and J.R. Heringa, J. Phys. A: 
Math. Gen. 28, 6289 (1995). 

E. Luijten, Interaction Range, Universality and Upper 
Critical Dimension (Delft University Press, Delft, 1997). 
E. Luijten and H.W.J. Blote, Phys. Rev. Lett. 76, 1557 

(1996). 

L.D. Landau and E.M. Lifshitz, Statistical Physics, (Perg- 
amon Press, London, 1959). 



11 



TABLES 



TABLE I: Variation of the critical temperatures and densities with - critical 

temperature obtained from the intersection of the cumulants Ql{P)', Tc{L) and pc{L) 
- fit of simulation data at 1^ = 8000 to the universal MF distribution; and pf - 
SCOZA (with a = 0.01); ri^^ and ' - MFL2. For statistical error see text. 



e 




Tc{L) 




rpS 

C 


pf 




fJc 


100.0 


113.5 


113.4 


0.250 


113.2 


0.249 


113.2 


0.249 


100/9 




12.62 


0.251 


12.60 


0.249 


12.60 


0.249 


1.0 




1.172 


0.242 


1.163 


0.243 


1.173 


0.243 


0.10 


0.1596 


0.1599 


0.190 


0.148 


0.176 


0.1604 


0.193 


0.05 


0.1054 


0.1062 


0.151 


0.0972 


0.0959 


0.1066 


0.135 


0.02 


0.0742 


0.0748 


0.111 


0.0825 


0.0225 


0.0850 


0.032 



TABLE IL Variation of the critical temperatures and densities with system size and 
Tc{L) and Pc{L) - fit of simulation data at volume V to the universal MF dis- 
tribution; pi{L) - from rule of rectilinear diameter (see text). For statistical error see text. 



e 


V 








0.020 


1000 


0.07595 


0.112 


0.102 


0.020 


4000 


0.07506 


0.111 


0.108 


0.020 


8000 


0.07477 


0.111 


0.108 


0.050 


1000 


0.10709 


0.153 


0.146 


0.050 


4000 


0.10620 


0.151 


0.147 


0.050 


8000 


0.10625 


0.151 


0.143 


0.100 


1000 


0.16125 


0.189 


0.183 


0.100 


4000 


0.16015 


0.189 


0.186 


0.100 


8000 


0.15993 


0.190 


0.187 


1.0 


1000 


1.181 


0.244 


0.244 


1.0 


4000 


1.174 


0.244 


0.243 


1.0 


8000 


1.172 


0.242 


0.239 


100.0 


1000 


114.0 


0.250 


0.249 


100.0 


4000 


113.5 


0.245 


0.248 


100.0 


8000 


113.4 


0.247 


0.246 



TABLE IIL Values of critical indices Q* , j/t, and yi obtained by the fit of Ql{P) [cf. Eq. 
MOIl l for difi'erent values of The fits have been performed with eight parameters. 



e 


Q* 


P? 


yr 


Vi 


0.02 


0.53±0.05 


13.47±0.1 


1.54±0.1 


-1.49±0.5 


0.05 


0.51±0.05 


9.48±0.1 


1.56±0.1 


-1.49±0.5 


0.10 


0.47±0.05 


6.27±0.1 


1.45±0.1 


-1.57±0.5 
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FIGURES 
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FIG. 1: Phase diagram in the (T, p)-plane of the system inves- 
tigated in the present study for = 1. Results were obtained 
via the improved MF theory (both one- and two-loop level) 
and SCOZA (for different values of a - see text); lines as 
indicated in the legend. In addition, a line connecting the 
respective critical points is drawn. 




logx 



FIG. 2: (Effective) critical exponent 7 as a function of r 
(for the definition see text) of the system investigated in the 
present study for = 0.16 as obtained from SCOZA for dif- 
ferent values of a, lined symbols as indicated in the legend. 
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FIG. 3: Variation of the critical temperature Tc with strength 
of the Kac potential: MF2L - squares and dashed line; 
SCOZA - triangles up and dashed line; solid circles - MC 
simulations for V = 1000 (see text). The solid line is a linear 
fit of the simulation data. Diamond and triangle down indi- 
cate the critical temperatures of the RPM ( ^^=0) obtained by 
the MSA approximation— and simulation—*—, respectively. 
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FIG. 4: Same as Fig. 3 but for the critical density pc- The 
MC results for the RPM are from Refi— (triangles up) and 
Reff^ (triangles right), respectively. 
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FIG. 5: Comparison between simulation and theoretical esti- 
mates for the coexistence curve at = 100/9 in the vicinity 
of the critical point. The symbols represent MC results for 

the different volumes: V = 1000 (circles), V — 4000 (squares), 
V = 8000 (triangles); MF2L (dashed line); SCOZA (fuU hue). 



FIG. 7: Coexistence curve at = 0.05 in the vicinity 
of the critical point. The symbols represent MC results 
for the different volumes: V = 1000 (circles), V = 4000 
(squares), V = 8000 (triangles); MF2L (dashed line). The 
inset shows SCOZA (full line) and the one-loop approxima- 
tion (diamonds). 
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FIG. 6: Coexistence curve at = 1 in the vicinity of the crit- 
ical point. The symbols represent MC results for the different 
volumes: V = 1000 (circles), V = 4000 (squares), V = 8000 
(triangles); MF2L (dashed line); SCOZA (full hne); one-loop 
approximation (dash-dotted line). 
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FIG. 8: Coexistence curve at ^'^ — 0.02 in the vicinity of 
the critical point. The symbols represent MC results for the 
different volumes: V = 1000 (circles), V — 4000 (squares), 
V = 8000 (triangles); MF2L (dashed hne); SCOZA (fuU line); 
one- loop approximation (diamonds). 
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are fits of the simulation data to Eq. 140^ . 
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FIG. 10: Fourth order cumulant Ql{P) as a function of in- 
verse temperature /3 at volumes V = 1000 (solid circles), 4000 
(squares), and 8000 (triangles) at — 0.02. The straight 
lines are fits of the simulation data to Eq. 1401 . 



